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Abstract 

In this paper wc introduce the concept of Lateral Trigger Probability (LTP) 
function, i.e., the probability for an extensive air shower (EAS) to trigger an in- 
dividual detector of a ground based array as a function of distance to the shower 
axis, taking into account energy, mass and direction of the primary cosmic ray. 
We apply this concept to the surface array of the Pierre Auger Observatory 
consisting of a 1.5 km spaced grid of about 1600 water Cherenkov stations. Us- 
ing Monte Carlo simulations of ultra-high energy showers the LTP functions are 
derived for energies in the range between 10^^ and 10^^ eV and zenith angles 
up to 65°. A parametrization combining a step function with an exponential 
is found to reproduce them very well in the considered range of energies and 
zenith angles. The LTP functions can also be obtained from data using events 
simultaneously observed by the fluorescence and the surface detector of the 
Pierre Auger Observatory (hybrid events). We validate the Monte-Carlo results 
showing how LTP functions from data are in good agreement with simulations. 

Keywords: Ultra-High Energy Cosmic Rays, Pierre Auger Observatory, 
Extensive Air Showers, Trigger performance. Surface detector. Hybrid 
detector. 
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1. Introduction 

The Pierre Auger Observatory has been eonceived to study the origin and 
the nature of uhra high-energy cosmie rays. Because of the scarcity of the flux 
at the highest energies, their direct measurement from space is technically un- 
feasible and the use of very large detectors is required at the ground. What can 
be observed is the extensive air shower of secondary particles produced in the 
propagation through the atmosphere. The Pierre Auger Observatory is located 
near Malargiie, Argentina, at 1400 m a.s.l. and it employes two independent and 
complementary measurement techniques The surface array (SD), consisting 
of about 1600 water Cherenkov detectors on a triangular grid of 1.5 km spacing 
covering an area of approximately 3000 km^ , records the secondary particles at 
the ground and thus samples their lateral density distribution. The fluorescence 
detector (FD), consisting of 24 telescopes at four sites, overlooks the surface 
array and observes the longitudinal profile of air showers by collecting the flu- 
orescence light emitted along the path through the atmosphere Q. Unlike the 
surface detector array with its nearly 100% duty cycle, the FD can only operate 
on clear and moonless nights giving an overall duty cycle of about 13% [3[. As 
a consequence, only a fraction of showers are observed by both detectors. For 
these events, called hereafter "hybrid", the combination of information from 
the surface array and the fluorescence telescopes enhances the reconstruction 
capability. Energy and direction reconstruction accuracy of hybrid events is in 
fact better than the one the SD and FD could achieve independently. 

One of the main goals of the Pierre Auger Observatory is to measure the flux 
of cosmic rays at the highest energies. This task relies on an accurate determi- 
nation of the detector exposure for SD-only [4j] and hybrid operation modes. 
The hybrid exposure is calculated using the simultaneous simulation of FD and 
SD response. Besides the dependence on energy and distance to an FD-site, 
the hybrid exposure is influenced by several factors including the atmospheric 
conditions, the trigger status of all active detectors and their instantaneous data 
taking configuration. The calculation of the SD response is based on the deep 
knowledge of the array capability to trigger once a shower with a given energy 
and zenith angle hits the ground. Since the trigger in an EAS array is always 
a combination of trigger states of neighboring detectors, the acceptance of any 
EAS array is directly connected to the probability that an individual detector 
triggers when a shower lands at a certain distance from it. This defines the con- 
cept of "Lateral Trigger Probability" function. This function has been used as a 
powerful tool for simulations in the analysis for the measurement of the hybrid 
energy spectrum and of the atmospheric depth at shower maximum 

The trigger design of the Auger surface detector is described in detail in [3| ■ 
Each water Cherenkov detector of the surface array has a 10 m^ water surface 
area and 1.2 m water depth, with three 9 in. photomultiplier tubes (PMTs) 
looking through optical coupling material into the water volume, which is con- 
tained in a Tyvek® reflective liner. The signals provided by each PMT are 
digitised by 40 MHz 10-bit Flash Analog to Digital Converters (FADCs) 
The achieved dynamical range is sufficient to cover with good precision both the 
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signals produced in the detectors near to the shower axis 1000 particles//is) 
and those produced far from the shower axis (~ 1 particle/^s). We recall here 
the basic structure of the used trigger algorithms. The two first levels (Tl and 
T2) are formed at each surface detector. Each trigger level can be divided in 
two modes, a threshold trigger (TH) and a time-over-threshold trigger (ToT). 
The first level threshold trigger (TH-Tl) requires the coincidence of the signals 
from the three PMTs equipping each station, each PMT signal being above 
1.75 "Vertical Equivalent Muon" (VEMf]. The TH-Tl trigger is used to reduce 
the rate due to atmospheric muons to ^100 Hz and can reach the second level, 
TH-T2, when the peak signal reaches at least 3.2 VEM in coincidence between 
3 PMTs signals, further reducing the rate to ^ 20 Hz. The second mode, the 
ToT, requires at least 13 time bins (i.e. more than 325 ns) in 120 FADC bins 
of a sliding window of 3 ^s to be above a threshold of 0.2 VEM in coincidence 
in 2 out of 3 PMTs. Time-over-threshold trigger stations are automatically 
promoted to the second level. The threshold trigger is especially efficient at 
detecting strong narrow signals, mostly encountered in horizontal showers or 
close to the axis of vertical showers. On the other hand, the ToT is intended to 
select sequences of small signals spread in time. This is typical of low energy 
vertical showers dominated by an electromagnetic component or of high energy 
showers triggering stations at large distance from the shower axis because of 
muons produced high in the atmosphere. 

Higher level triggers are obtained by requiring the spatial and temporal 
coincidence of at least three stations satisfying the T2 conditions. In particular, 
for zenith angles below 60°, the full efficiency for SD is reached at 10^^'^ eV 
In addition, if at least one FD telescope triggers in coincidence with one second 
level trigger station, a hybrid trigger is formed. Since FD has a lower energy 
threshold, hybrid events are also detected below the minimum energy for an 
independent SD trigger. For zenith angles below 60°, the hybrid detector reaches 
nearly full efficiency at 10^® eV 

In section 2, the concept of a Lateral Trigger Probability (LTP) function 
is formalized and applied to the particular case of the surface detector of the 
Pierre Auger Observatory. In section 3, the LTP functions for a single time- 
over-threshold trigger station are derived and parametrized for different primary 
particles (proton, iron, photon) and their dependence on energy and zenith angle 
is explored for zenith angles up to 65° and for energies between 10^^ and 10^^ eV. 
This energy range is relevant as it covers the interval in which the SD-only and 
the hybrid detection mode become fully efficient. The dependence on the choice 
of the hadronic interaction models is also discussed in section 3. In section 
4, hybrid data arc finally used to validate the simulation and to estimate the 
impact of weather conditions on the observed efficiency. The LTP functions are 



^ The distribution of measured light due to atmospheric muons produces a peak in the 
PMT charge distribution, Qy'em i° short), as well as a peak in that of the pulse 

height, ^vEM' both of them being proportional to those produced by a vertical through-going 
muon [3]. 
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found to reproduce very well the detector response over a wide range of energies 
and zenith angles. 

2. Concept of Lateral Trigger Probability 

The trigger probability of a single water Chercnkov detector depends on 
several independent physical parameters : i) the characteristics of the primary 
cosmic ray that initiates an air shower, e.g., its energy and mass, ii) the type 
and geometry of the detector used to observe air showers (in the following we 
will only study water Chercnkov detectors used for the surface detector array of 
the Pierre Auger Observatory), iii) the trigger condition used to detect a signal 
from air showers iv) the geometry of the incoming shower, e.g. its incidence 
zenith angle and position with respect to the detector. To formalize these de- 
pendencies we define the Lateral Trigger Probability function IS.A,E.e,Tr{'r, 0) as 
the probability to trigger on an air shower induced by a primary particle of en- 
ergy mass A and zenith angle 0. Here, r and (j) £^re the radial coordinates 
of the single detector in the plane normal to the shower axis (shower frame). 
Using a trigger condition Tr, this probability is simply given by: 

Ni 

^Am{r.<t>) = j^-^^- (1) 

where A'^i and A'o arc respectively the number of triggered and un-triggered 
detectors with coordinates r and (j) in the shower frame. 



3. Simulations 

The LTP functions have been derived using detailed simulations of the EAS 
development and of the detector response. The simulation sample consists of 
about 15000 CORSIKA Q showers (proton, iron and photon primaries) with 
zenith angle distributed as sin cos 6^ {0 <65°) and energies ranging between 
IQ-*^^ and 10^^ eV in steps of 0.25 in the logarithmic scale. A "thin sam- 
pling" mechanism at the level of 10~^ (optimal thinning) is applied following 
the standard method used for CORSIKA simulation with energies larger than 
10^6 eV 0. The showers have been generated with the models QGSJETII ^ 
and FLUKA [llj for high and low energy hadronic interactions. 

In the simulation, the position of the shower core (i.e. the intersection of the 
shower axis with the ground) is uniformly distributed over the surface array and 
each shower is used 5 times, each time with a different core position, in order 
to cover different areas of the array and explore all the detector configurations. 
The surface detector response is simulated using GEANT4 [l^l and adopting the 
sampling procedure to regenerate particles in a ground detector from thinned 
air shower simulations as described in 1131. The entire detector simulation is 



carried out within the framework provided by the Auger Offline software [14 1 . 

The trigger status of SD stations is inspected within a radius of 3 km from 
the shower axis and the Lateral Trigger Probability is then derived according to 
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eq. [TJ At distances larger than 3 km, the trigger efficiency is negligibly small for 
the class of events studied in this paper. All trigger modes of the surface detec- 
tor are simulated in detail at all levels. However, for events with zenith angles 
below 65°, the majority of the stations forming a second level trigger satisfy 
the ToT condition. In particular, for the considered zenith angles, the fraction 
of TH-T2 trigger stations not being also ToT is about 1%, approximately inde- 
pendent of the energy. Thus, we focus the analysis on the ToT stations. 

The lateral trigger probability for a ToT station is shown in Fig. [T] at a given 
energy and for different ranges of the cosine of the zenith angle 9. The max- 
imum effective distance for detection increases with energy and, for a given 
energy, with the cosine of the zenith angle, i.e. events with larger zenith angle 
tend to trigger less due to the attenuation of their electromagnetic component. 
For moderately inclined showers, an asymmetry is expected in the signal de- 
tected in the stations placed at the same distance to the shower axis but with 
different azimuth in the shower frame [lH . Indeed, secondary particles arriving 
earlier traverse less atmosphere and are less attenuated than the late ones. As a 
consequence, early stations may exhibit larger trigger probabilities and produce 
larger signals. Actually, for zenith angles below 65°, this effect has been found 
to have a quite low influence on the trigger probability, only noticeable above 
30° (in simulations as well as in the data). In the following we consider LTP 
functions averaged over all azimuths in the showers frame. A more detailed 
treatment including the azimuthal dependence does not introduce measurable 
differences for acceptance calculations. 

A fit combining a step function (close to the axis) with an exponential (fur- 
ther away) reproduces reasonably well the full simulated data set. The form of 
the fit function used is: 



where Rq, AR and C are free fit parameters, with Rq being the distance where 
LTP is equal to 0.5. A fit performed according to eq. [2]is superimposed on each 
plot shown in Fig. [T] As an example, the ToT trigger probability at energy 
£^=10^^ eV and for two angular bins (vertical showers on the left and showers 
with larger zenith angle on the right) is shown in Fig. [5] the exponential can 
reproduce very well the tail of the probability distribution at large distances 
from axis, in particular for inclined events. The dependences of fit parameters 
Rq, AR and C on energy and zenith angle can be parametrized by quadratic 
polynomials in the variables cos 9 and \ogiQ{E/eV). The corresponding coeffi- 
cients are tabulated in the Appendix for proton, iron and photon primaries. In 
Fig. [3l the ToT trigger probability from parametrization has been superimposed 
on the simulation (proton primary, all zenith angles up to 65° are merged). The 
comparison is performed as in the following. For each simulated event, i.e. for 
a certain primary, energy and arrival direction, the LTP is calculated using the 
parametrization (lines) and shown together with the full simulation (points). 
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Figure 1: Lateral Trigger Probability from simulations (proton primary) for a 
ToT station at a given energy, from 10^^ eV up to 10^^ eV in steps of 0.5 in 
the logarithmic scale. Different bins of cos 6 are also shown together with a fit 
performed according to eq. |2l superimposed as a continuous line. 




Figure 2: Fit made with a step function in proximity of the shower axis (con- 
tinuous line) and by an exponential at larger distances (dashed line). The ToT 
probability is shown for vertical (left) and inclined (right) showers at energy of 
10^^ cV. 
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The agreement is remarkably good in the entire energy range for proton (shown 
in the figure) and for iron and photon primaries. 

3.1. Dependence on primary mass 

The detector response to showers induced by different primary particles is 
shown in Fig. 01 for two classes of events, vertical (0° < 6 < 38°) on the top 
and moderately inclined (38° < 6 < 65°) on the bottom. Because of their 
larger number of muons, showers induced by iron nuclei provide a higher trigger 
capability at larger distances than those induced by protons, for all zenith angles. 
However, the difference between proton and iron is too small to give any hint for 
mass composition analysis. On the other hand, the LTP functions for photon 
primaries differ sensibly from those of hadrons (they vanish at shorter distances, 
about 500 m less at an energy of 10^^ eV). This is a consequence of the structure 
of the lateral distribution of photon showers, i.e. at a given energy, their effective 
footprint at the ground is smaller than the one of hadrons. Moreover, in photon 
showers there is a much smaller number of muons. 

It is worth noting that the energy threshold corresponding to full efficiency 
for SD, derived from data and simulation in ref. has been found to be 
compatible with the expectation for hadronic primaries. 

3.2. Dependence on hadronic interaction model 

Different choices of high energy interaction models influence the simulation of 
shower development and could affect the expected trigger efficiency. The depen- 
dence of the Lateral Trigger Probability on the assumptions for the hadronic 
interaction model has been investigated using a sample of simulated showers 
(proton) produced with SIBYLL [l6| . As shown in Fig. [SJ the LTP functions 
derived with the two hadronic interaction models differ only at large distance 
from the shower axis, in a range where the efficiency degrades rapidly. In this 
region, SIBYLL gives a lower LTP since this model predicts on average a smaller 
number of muons. Those differences are however too small to imply an observ- 
able impact on the detector acceptance. 

4. LTP functions from data and comparison with simulation 

The LTP functions can be derived from data by calculating the ratio of trig- 
gered to active stations within a given distance from the reconstructed shower 
axis. While doing this, the actual surface detector configuration must be accu- 
rately taken into account as a function of time. In addition, only high quality 
data are selected to avoid biases due to mis-reconstructed energies and/or ar- 
rival directions. The use of hybrid events allows to derive LTP functions also 
for energies below the threshold of an independent SD trigger. This is a benefit 
of the hybrid design that aims to fully exploit the distinctive potential offered 
by the Pierre Auger Observatory. Two years of hybrid data collected between 
June 2006 and May 2008 were used for this study. The events are selected as 
described in and this ensures an angular resolution of about 0.6° and a core 
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0° < zenith < 65° 




Figure 3: Lateral Trigger Probability for a ToT station as a function of station 
distance to shower axis and for different energies (proton primary) . The outcome 
of the parametrization is superimposed as a line. All zenith angles up to 65° 
are merged. 
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Figure 4: Lateral Trigger Probability for a ToT station. Proton, iron and 
photon primaries of energy 10"'^^ eV for two zenith angle ranges, 0° - 38° (top) 
and 38° - 65° (bottom). The outcome of the parametrization is superimposed 
as a continuous line. 
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Figure 5: Lateral Trigger Probability for a ToT station (zenith angle between 0° 
and 65°). Proton primary at energy of lO^^ eV with QGSJETII and SIBYLL. 
The ratio QGSJETII/SIBYLL is shown in the inset. 
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position determination better than 70 m. Further requirements on the goodness 
of the reconstructed longitudinal profile provide an energy resolution of about 
10% above 10-^^ eV and less than 15% at lower energies Q. 

The LTP measured from data is shown in Fig. [5] for different energy intervals. 
To verify the performance of the parametrization described in section [31 for 
each selected event, the LTP of any active station within 3 km from the shower 
axis is calculated using the reconstructed energy and direction. The predicted 
probability (dashed line) is then superimposed on data (points), see Fig. [Sj In 
this way, data are compared to simulation taking into account the actual status 
of the detector. The shaded area gives the interval of expected values assuming 
that data are pure proton (lower edge) or pure iron (upper edge) . A 50% proton 
and 50% iron mixed composition has been assumed for the parametrization 
(dashed line). The agreement is good over the entire energy range. This feature 
actually starts at very low energies, even below the range of full efficiency for 
the hybrid detection [3[. In this case, whereas in data only events with at least 
one SD ToT station are selected, in simulation also the events that did not 
trigger at all are taken into account in the calculation of the probability. As 
a consequence, the comparison between data and simulation could be biased. 
However, the good level of agreement actually reached reflects the fact that the 
hybrid detection is very close to fully efficient and the energy reconstruction 
remains reasonably good within the scope of this analysis down to energy of 
about IQi^-^ eV. 

For each energy interval considered, the agreement between data and simu- 
lation has also proven to hold in two zenith angle bands (0° - 38° and 38° - 65°), 
see Fig. [71 



Impact of weather effect on LTP 

The effect of atmospheric variations (in pressure, temperature and air den- 
sity) on extensive air showers development has been extensively studied with 



the surface detector data [17[. A significant modulation of the rate of events 



with the atmospheric variables, both on a seasonal scale 10%) and on a 
shorter time scale (^ 2% on average during a day) has been observed. This 
modulation is mainly explained as due to the change with the air density of the 
Moliere radius near ground thus influencing the trigger probability and the rate 
of events above a fixed energy. Hybrid data in the energy range around 10^^ eV 
have been used to investigate this effect on LTP. Data have been separated by 
season and are shown, together with the parametrization, for austral winter and 
austral summer, see Fig. [8l top panel. The ratio of summer and winter rela- 
tive to the parametrization is shown in the bottom panel. Results qualitatively 
match the expectation. Higher temperature at the ground, as for the austral 
summer, induces a reduction of the air density weakly enhancing the trigger 
probability at a given distance relative to all other seasons. Nevertheless the 
effect is almost negligible on the scale of the measurable trigger efficiency. 
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Figure 6: Comparison of simulation with hybrid data collected in two years. All 
zenith angles up to 65° merged. The energy intervals are 10^^'^ < E < 10^^'^ eV, 
IQi^-^ <E < 1018-2 eV, 1018-2 <E < lO^^-^ eV, IQi^-^ <E < lO^^-^ eV. 



5. Summary and Conclusions 

In the previous sections we have introduced the concept of Lateral Trigger 
Probability function as a tool to characterize the single detector trigger effi- 
ciency. We have derived LTP functions for the particular case of the surface 
detector of the Pierre Auger Observatory using simulations. We discussed their 
evolution with different physical parameters of air showers such as the energy, 
zenith angle and nature of the primary particle. We also investigated the impact 
of choosing different hadronic interaction models in the simulations. Further- 
more, we estimated the LTP functions at different energies and zenith angles 
using hybrid data and showed that seasonal effects are visible in the trigger 



probabilities retrieved from data as expected from previous studies [17|. 

The good agreement between simulations and data over a wide energy range 
(between 10^''-^ eV and lO^^ eV) demonstrates the accuracy of the different 
aspects of the simulation procedure (i.e. air shower, detectors and trigger sim- 
ulation) as well as the quality of the reconstruction obtained for hybrid data. 
These comparisons support and validate the use of simulated LTP functions 
in the estimate of the hybrid aperture described in Monitoring the LTP 
functions over a longer period of time can be used to study the long-term per- 
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Figure 7: Comparison of simulation with hybrid data collected in two years. 
Zenith angles are split in two ranges 0° - 38° (left) and 38° - 65° (right). 
From top to bottom the energy intervals are lO^'''^ < E < 10^^'^ eV, 
1017-^ <E < 1018-2 eV, 1018-2 <E < 1018-7 eV, IQi^-^ <E < IQi^-^ eV. 
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formance of the SD trigger for individual stations both above and below the 
acceptance saturation energy. 

As a final consideration, LTP functions can be derived at higher energies 
using SD-only data because, at energy above ^ 10^®-^ eV, despite the statis- 
tics of hybrids becoming small, the surface detector is fully efhcient and the 
geometrical reconstruction is accurate. As mentioned in the Introduction, the 
probability of a high level trigger for the surface detector is a combination of 
single detector probabilities. Hence LTP functions provide a robust and simple 
method to estimate the energy or zenith angle dependence of SD acceptance for 
any arbitrary configuration. This makes this technique a valuable tool to design 
other experiments and future enhancements of the Pierre Auger Observatory. 
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Appendix: LTP parametrization 

The LTP is fitted, as discussed in section |31 to the fohowing function: 

LTP(r)J ^-"^ (A.1) 

with Rq being the distance where the LTP is equal to 0.5. 
The dependences of fit parameters Rq, AR and C on energy and zenith an- 
gle can be parametrized by quadratic polynomials in the variables cos 9 and 
logio{E/eV). The corresponding coefficients are given for proton, iron and 
photon primaries {0° < 9 < 65°), separately. Concerning the accuracy of the 
parameters, a change at the level of (1 -^5)% propagates approximately linearly 
in the returned value of the parametrization. 



Proton showers 

The overall parametrization for proton primaries (0° < 9 < 65*^ 
rized in the following matrix equation: 
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Iron showers 

The overall parametrization for iron primaries (0° < 9 < 65°) is summarized 
in the following matrix equation: 
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Photon showers 

The overall parametrization for photon primaries (0° < 6 < 65°) is summa- 
rized in the following matrix equation: 
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